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ABSTRACT 

Using new photometric and spectroscopic data in the fields of nine strong gravitational lenses that 
lie in galaxy groups, we analyze the effects of both the local group environment and line-of-sight 
galaxies on the lens potential. We use Monte Carlo simulations to derive the shear directly from 
measurements of the complex lens environment, providing the first detailed independent check of the 
shear obtained from lens modeling. We account for possible tidal stripping of the group galaxies by 
varying the fraction of total mass apportioned between the group dark matter halo and individual 
group galaxies. The environment produces an average shear of 7 = 0.08 (ranging from 0.02 to 0.17), 
significant enough to affect quantities derived from lens observables. However, the direction and 
magnitude of the shears do not match those obtained from lens modeling in three of the six 4-image 
systems in our sample (B1422, RXJ1131, and WFI2033). The source of this disagreement is not clear, 
implying that the assumptions inherent in both the environment and lens model approaches must be 
reconsidered. If only the local group environment of the lens is included, the average shear is 7 = 0.05 
(ranging from 0.01 to 0.14), indicating that line-of-sight contributions to the lens potential are not 
negligible. We isolate the effects of various theoretical and observational uncertainties on our results. 
Of those uncertainties, the scatter in the Faber- Jackson relation and error in the group centroid 
position dominate. Future surveys of lens environments should prioritize spectroscopic sampling of 
both the local lens environment and objects along the line of sight, particularly those bright (I < 
21.5) galaxies projected within 5' of the lens. 
Subject headings: gravitational lensing: strong 



1. INTRODUCTION 

Analyses of strong gravitational lenses have been use- 
ful in probing cosmological parameters such as Hp 
(e.g.|Refsdallll964tlKeeton fc Kochaneklll997l;ISaha et alJ 
12006c lOguril |2007|) . constraining p roperties of the 
dark matter halos of gal axies (e.g. iKoopmans et al.1 
l20"06t IBarnabe et all l200l . and uncovering su b struc- 
ture in those halos Ce.g. IMao fc Schneider! 119981 : 
iMetcalf fc MadavfeOOlHDalal fc Kochaneldl2002¥ . These 
studies require accurate models of the lens potential, 
which can have contributions not only from the main 
lens galaxy, but from other objects at the lens redshift 
or along the line of sight. Indeed, lens models often 
req uire environmental terms, representing a tidal shear 
H: IKeeton et al.lll997h and perhaps higher-order effects 
(lKochaneld(l99Tl : IKeeton fc Zabludofi1[200l IPadelv et al.1 
12010ft . inorder to yield a good fit to the observed im- 
age positions and flux ratios. If we could measure lens 
environments, we would have independent and direct de- 
terminations of the shear to compare with lens models 
as a test of their results. 

Past studies have suggested that the local lens en- 
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vironment can have a non-negligible effect on the lens 
model-de rived shear of a sys tem. Using galaxy demo- 
graphics, Kee ton et al.1 (|2000T ) estimate that at least 25% 
of strong lenses are in group or cluster environments 
that could cause strong perturbati ons in the lens po- 
tentia l. With N-body simulation s, iHolder fc Schechterl 
(2003) and lDalal fc Watsonl (|2004l ) compute the environ- 
mental shear at the likely positions of lens galaxies, al- 
though they reach somewhat differe nt conclusions about 
how s trong the shear should be. [Holder fc Schechterl 
(2003) find an expected value of 7 ~ 0.11, which is sim- 
ilar to t he amount of shear need ed in some lens models, 
whereas iDalal fc Watsonl (|2004f ) find a lower expected 
value of 7 w 0.03. While useful, such statistical stud- 
ies do not capture the richness and possible diversity of 
individual lens environments. 

Direct calculations of the effects of lens environ- 
ments requir e extensive observation a l data. Sev- 
eral st udies dFassnacht fc Lubinl l2002t iFassnacht et al.1 



2006t iMomcheva et all 120061: [Williams et al.1 1200a 
Auger et all 120071: iMoustakas et all 120071: lAuger et al l 



2008t iFaure et al l 120091 120101 iMcKean et all I2010T 



Suvu et aLl I2010T ) have used observational data to es 



timate the effects of environment on the potentials of 
strong lenses. However, these studies did not calculate 
the shears induced by the environment to a level where 
they could draw conclusions about the most significant 
sources of the lens potential perturbation, nor were they 
always able to compare the calculated shears to the re- 
sults of lens models. Using spectroscopic and photomet- 
ric data (Momcheva et al. in preparation; Williams et al. 
in preparation) for nine strong gravitational lenses that 
lie in groups, we measure the environments both around 
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the lens and projected along the line of sight (LOS) to 
the lens to determine the extent to which they affect the 
lensing potential. We refer to the combined local environ- 
ment plus LOS perturbers as the full lens environment. 
Agreement between the shears inferred from the environ- 
ment and those derived from lens models would demon- 
strate that the environment is responsible for the pre- 
viously unexplained large shears found for some lenses. 
Any disagreements would point to problems in either the 
environment treatment or lens models. 

Objects projected close to the lens can contribute to 
the lens potential, even if they are not physically associ- 
ated with the lens. The magnitude of the LOS perturba- 
tions decrease with both increasing separation from the 
lens (radial dis tance and redshift) and decreasing mass of 
the perturber (Mo mcheva et al.1l2006| ). We must charac- 
terize the effects of the LOS environment as a function of 
both projected separation and apparent magnitude in or- 
der to estimate the point at which LOS perturbers' con- 
tributions to the lens potential become negligible. This 
will better inform the observational strategies of future 
surveys of lens environments by placing estimates on the 
size and depth of a spectroscopic survey that will suffi- 
ciently characterize the most significant contributions to 
the lens perturbations. 

Shears derived from measurements of the lens environ- 
ments may be affected by a number of systematic and 
random uncertainties in the theoretical formalism behind 
our methodology, as well as by observational errors and 
incompletenesses in our data. The observational uncer- 
tainties for which we explicitly account include uncer- 
tainty in the Faber- Jackson relation, errors in the pro- 
jected group centroid position and velocity dispersion, 
and magnitude errors in our photometry. The theoretical 
uncertainties include the apportionment of mass between 
the group dark matter halo and the halos of individual 
group galaxies, the form of the density profile for the 
group dark matter halo, and scatter in the concentration 
parameter of the group halo. By quantifying the relative 
importance of these effects, we can determine those for 
which we must account in future lens surveys, as well 
as improve the theoretical assumptions used to calculate 
the lens environments. 

This paper is organized as follows. Our sample of lens 
systems, along with details of our spectroscopic and pho- 
tometric catalog, is described in § [2j The lens modeling 
of the 4-image lenses is in § [3] In §|4j we discuss our meth- 
ods for measuring the lens environments, including our 
handling of incompleteness in our spectroscopic catalogs 
and additional sources of error introduced by our theo- 
retical assumptions and observational uncertainties. We 
also describe our methodology for quantifying environ- 
ment effects on the lens potential. In § [5j we describe our 
treatment of the local group environment and perturbers 
along the line of sight to the lens, as well as the shear 
due to the full lens environment, including both local 
group and LOS perturbers simultaneously. We present 
our main results in § [6] and summarize our conclusions 
in § In the appendices, we present the details of the 
shear arising from assumptions about the group and indi- 
vidual galaxy halo mass distributions (Appendix [A|, our 
calculation of the truncation radii of group galaxies (Ap- 
pendix [B]) , and the effects of radial and luminosity cuts 
on line-of-sight shear (Appendix [C]). Throughout this 



paper, we assume a flat cosmology with Q m = 0.274, fif, 
= 0.045, n A = 0.726, and H = 71 km s" 1 Mpc" 1 . 

2. THE SAMPLE 

Our sample consists of nine strong gravitational 
lenses chosen from the full sample of 26 lenses 
analyzed by Momcheva et al. (in preparation). 
These nine systems are those that were deter- 
mined to be in galaxy groups: Q ER 0047- 
2808 (hereafter Q0047; iWarren et aTT 119961 11998ft. 
HE 0435- 1 223 (hereafte r HE0 4 35; iWisotzki et atl 



120001 [20M IMorgan et ail ifM IQfek et al.l I27)06ft~ 



MG J0751+2716 (hereafter MG 0751; iLehar et all 
[1991 [19971 iTonrv fc Kochanekl fl999ft . P G 1115+080 
/hereafter PG11 15: iWevmann et alj |1980t Kundic et all 
Il997al ITonrv! Il998ft RX J1131-1231 (h ereafter 
RXJ1131~ lSluse et al.l [20031 IMorgan etalll2006h. HST 
J14113+5211 (herea fter HST14113; iFischer et all 
1998HLubin et al.ll2000D, B1422+231 (h e reafter B1422; 



Patnaik et al.l Il992t llmpev et all 119961 : iKundic et all 



1997b:" TKing et all 119991), MG J1654+1346 (hereafter 



MG1 654~ TLangston et all 119881 Il989t iKochanek et al.l 
1200(1 . and WFI J2033-4723 (hereafte r WFI2033; 
Morgan et all [200l lEigenbrod et~aTI [2001 IQfek et alJ 
2006|) B 

Spectroscopic observations of four of these lenses 
and th eir environments are detailed in iMomcheva et all 
(2006), with further spectroscopic data on all nine sys- 
tems in Momcheva et al. (in preparation). We also have 
two-band phot o metri c catalogs of the lens fields from 
iWilliams et al.l (|2006ft and Williams et al. (in prepara- 
tion). Figures Q] and [2] show the projected spatial distri- 
bution of objects in our catalogs. The lens properties are 
in Table [TJ and the host group properties are in Table [5J 

The spectroscopic data were taken over multi- 
ple observing runs with the Low Dispersion Survey 
Spectrograph-2 (LDSS-2; lAllington-Smith et al.l 119901 1 
and LDSS-3 on the 6.5 m Magellan 2 (Clay) telescope, 
as well as with the Inam ori Magellan Aerial Cam era 
and Spectrograph fIMACS; lBigelow fc Dresslerll2003T ) on 
the 6.5 m Magellan 1 (Baade) telescope, both at Las 
Campanas Observatory. Additional spectroscopic data 
were taken with the H ectospec multi-object spectrograph 
(jFabricant et al.ll2005ft on the 6.5 m Multiple Mirror Tele- 
scope on Mt. Hopkins. The spectroscopic target selec- 
tion prioritized objects brighter than / = 21.5 and within 
5' of the lens in projected separation. We plot magnitude 
histograms for objects within 5' of each lens galaxy in 
Figure [3 with our spectroscopic limit of / = 21.5 shown 
for comparison. 

6 T hree additional le ns systems, SB S1520+530 HBurud et al.l 
[2001: lAuger et alj 12008ft. B1600+434 Uaunsen &: Hiorthl I1997F 
Fassnacht fc Cohenl IT998: Koopma ns et al .l 119981 ; lAuger et al.l 
120071 ), and B2114+022 lIKing et al.l 119991 ; lAugusto et al.l I200JT 

may also lie in galaxy groups. However, we exclude these from 
our analysis. SBS1520 has an uncertain lens redshift, which affects 
its membership in a potential host group. The published veloc- 
ity dispersion of the group thought to be associated with B1600 
is less than 100 km s _1 and determined from only five member 
galaxies (Momcheva et al. in preparation), which makes it very 
uncertain. B2114 lacks a measured source redshift a nd is being 
stron gly lensed by two galaxies at different redshifts ( Ch ae et al.l 
2001), complicating the analysis of its lens potential. Furthermore, 
all three of these systems are two-image lenses, which provide little 
if any constraint on the shear. Therefore, these systems are lower 
priority than those for which we do have lens models. 
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Aa (arcmin) 

Fig. 1. — Projected spatial distribution of galaxies in the lens fields. North is up and East is to the left. The fields are centered on the lens 
(purple star). The black cross represents the projected group centroid and its error. Also shown are the spectroscopically confirmed group 
members (open blue circles), spectroscopically confirmed line-of-sight (LOS) objects (green squares), photometric red sequence galaxies 
down to / = 21.5 (red circles), and the remaining photometric galaxies down to I = 21.5 (black dots). The radii of the circles representing 
the group members are scaled by their Einstein radii, Rg. The lens galaxy is scaled in a similar manner. The size of the black circle in 
the upper right corner represents a galaxy with an internal velocity dispersion of a = 200 km s . The area of the points representing the 
lens and group members scales as luminosity within a panel, but is not consistent between panels because of the D^g/Dg scaling between 
a and Re- The bar in the lower left corner corresponds to 1 Mpc at the lens redshift. 



The ph otometric data are fr om imaging described in 
detail in iWilliams et al.l (|2006f ) and Williams et al. (in 
preparation). Wide- field imaging is obtained for each 
field with the Mosaic imagers on the Kitt Peak and CTIO 
4-meter telescopes in two photometric bands, Cousins / 
and either Johnson V or Cousins R. Total exposure times 
vary in order to ensure high completeness for / < 21.5. 
We obtain at least one image of each field on a photomet- 
ric night. Exposures obtained on non-photometric nights 
are corrected to the standard system using local stellar 
calibrators derived from the photometric images. Object 
detection and phot ometric measurements are performed 
using SExtractor ()Bertin fc Arnoutill996() . 

/-band magnitudes are MAG_AUTO output from SEx- 



tr actor, which are calculated in a manner similar to that 
of iKronl |l980). Total magnitudes, which we later use 
to calculate gala xy velocity distributions from the Faber- 
Jackson relation (jFaber k, Jackson! 1 9 76t l , are determined 
as follows. For each field, well-isolated galaxies with no 
SExtractor photometric flags are selected. We measure 
magnitudes in very large apertures, determine their off- 
set from the MAG_AUTO values, and use the mean off- 
set to correct all objects in the catalog. These offsets are 
small (~ 0.05 mag) with small dispersions (~ 0.03 mag). 
All magnitudes used in this paper are Kron magnitudes 
unless indicated otherwise. 

SExtractor fails to detect some galaxies, particularly 
those in the halos of relatively bright stars. We manually 
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Fig. 2. — Same as Figure[T] but with a smaller field of view. The bar in the lower left corner of each panel corresponds to 100 kpc at the 
lens redshift. 



TABLE 1 
Gravitational Lens Properties 



Lens 


a (J2000) 


8 (J2000) 


# images 




^source 


Re(") 


Q0047 


12.425 


-27.874 


ring 


0.484 


3.60 


1.35 


HE0435 b 


69.562 


-12.287 


4 


0.455 


1.69 


1.21 


MG0751 


117.923 


27.276 


ring 


0.350 


3.20 


0.35 


PG1115 b 


169.571 


7.766 


4 


0.310 


1.72 


1.16 


RXJ1131 b 


172.965 


-12.533 


4 


0.295 


0.66 


1.90 


HST14113 


212.832 


52.192 


4 


0.464 


2.81 


0.90 


B1422 b 


216.159 


22.934 


4 


0.337 


3.62 


0.84 


MG1654 


253.674 


13.773 


ring 


0.253 


1.74 


1.05 


WFI2033 b 


308.425 


-47.395 


4 


0.661 


1.66 


1.16 



Note. — Lens redshifts and image separations from CASTLeS and Rusin et al.] 120031 ). 



^Re estimated to be half image separation of lens. 
b Lens has measured time delay. 
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TABLE 2 
Host Group Properties 



Lens 


# members a 


a cen (J2000) 


5 cen (J2000) 


acen error (' 


) b Seen error (' 


') b Group oy (km s" 1 ) c 


Group / RS 


Q0047 


9 


12.426 


-27.862 


0.26 


0.48 


348t^ 




HE0435 


11 


69.562 


-12.290 


0.50 


0.66 


p- 99 +106 




MG0751 


29 


117.908 


27.289 


0.50 


0.50 


518± 8 f 


5/11 


PG1115 


13 


169.568 


7.765 


0.35 


0.26 


390 + |° 


3/11 


RXJ1131 


27 


172.896 


-12.570 


0.64 


0.57 


429±J| 8 


7/13 


HST14113 


41 


212.818 


52.204 


0.43 


0.46 


656±gg 


17 / 31 


B1422 


17 


216.176 


22.931 


0.42 


0.39 




6/14 


MG1654 


8 


253.663 


13.791 


0.30 


0.38 




2/11 


WFI2033 


14 


308.449 


-47.365 


0.47 


0.33 


498l 8 7 | 





Note. — Tabulated values calculated using methods in Momchcva ct al. (in preparation), consistent with methods described in text. 
a Includes only spectroscopically confirmed members. 
b Errors calculated from bootstrap resampling. 

c Errors calculated from jackknifc resampling and bi-weight estimators as in Momchcva ct al. (in preparation). 

d (# spectroscopically confirmed group members on the red sequence) / (# red sequence galaxies with spectroscopy). Groups with no 
photometric red sequence at lens redshift are omitted. 
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Fig. 3. — Histogram as a function of apparent /-band magnitude 
for galaxies in the photometric sample (open) and spectroscopic 
subsample (hatched). Only objects projected within 5' of the lens 
are included. Our spectroscopic completeness drops off at I = 21.5 
(dashed li ne), t he limit at which we cut our sample for analysis 
(Appendix l02t . 

determine positions and aperture magnitudes for those 
objects with / < 21.5 and within 4' of the lens galaxy. 
Total magnitudes for these galaxies are determined in 
a similar fashion as described above. These corrections 
tend to be large (> 0.5 mag) and the scatter large (~ 0.3 
mag), as small apertures are used in order to minimize 
the significant shot noise from the bright star halos. 

We obtain the lens galaxies' image separations from the 
CfA-Arizona Space Telescope Lens Su rvey (CASTLeS) 
database and from iRusin et all (|2003ft . The HST data 
from CASTLeS are needed to separate the positions and 
fluxes of the lensed images from the lens galaxy because 
our photometry is not always able to resolve the individ- 
ual components. 

We exclude serendipitously observed stars from the 
spectroscopic catalog (z < 0.01). We also exclude high- 
rcdshift QSOs and AGNs (z > 1.5). For similar rea- 
sons, we exclude objects whose inferred /-band absolute 
magnitudes are brighter than —25 as galaxies this bright 
should be exceedingly rare (jBlanton et al.l [2001) given 
the volume of our survey. Of all the galaxies projected 
within 5' of a lens galaxy and brighter than / = 21.5 
across all nine fields in our spectroscopic sample, fewer 
than ~5% are removed by these cuts. There is one ob- 
ject in the field of B1422 at a=216.200, 5=22.9475 that 
is spectroscopically confirmed to be a group member but 
does not have robust photometry because it lies under 
a bleed trail. We use it in calculating the group cen- 
troid position and group velocity dispersion, but do not 
include it in our shear analysis. Given its projected dis- 
tance from the lens, it is unlikely to significantly affect 
the shear calculation. 

3. LENS MODEL SHEAR 

We model the six 4-image lenses in our sample to cal- 
culate the shears needed to fit the lens data. We consider 
as constraints the positions and fluxes of the lensed im- 
ages, along with the position of the lens galaxy. For 



HE0435, RXJ1131, HST14113, and WFI2033, we use 
position and relative flux data from CASTLeS, using 
the longest wavelength band available (typically WFPC2 
F814W or NICMOS F160W) to limit the effects of mi- 
crole nsing and variability . For B1422, we use radio data 
from lPatnaik et ail (|1999fl to get the relative image posi- 
tions and fluxes, and CASTLeS data to get the position 
of the lens galaxy relative to the images. For PG1115, 
we take the image posi tions from CASTLe S, but use the 
mid-IR flux ratios from Chi ba et al.l (|2005D . The position 
uncertainties are 3 mas except for the following cases: 8 
mas for the RXJ1131 lens galaxy, (7,38,9,9) mas for the 
lensed images (A,B,C,D) in HST14113, and 9 mas for 
the WFI2033 lens galaxy. We broaden the error bars 
on the fluxes to 10% to account for microlensing and/or 
variability, except for P G1115 where we us e the measure- 
ment errors reported by Chiba ct al. (2005), since mid-IR 
flux ratios should not be susceptible to those systematic 
effects. 

The models assume an ellipsoidal lens galaxy with a 
power law mass density profile, plus external shear. The 
power law index, a, is defined such that the enclosed 
mass projected within R is M{R) oc R a . We do not place 
any explicit constraints on the lens galaxy's Einstein ra- 
dius, ellipticity, and position angle, or on the amplitude 
and direction of the shear Q While the ellipticity and 
position angle have been measured for some lens galax- 
ies, we do not want to assume that the mass necessarily 
traces the light. We run Monte Carlo Markov Chains 
(MCMC) to sample the full range of allowed models. For 
the purpose of this paper, the key model output is the 
distribution of shear values obtained after marginalizing 
over all of the lens model parameters. We do not in- 
clude any external convergence in the lens models, since 
that cann ot be constrained due to the mass sheet de- 
gcneracy (|Falco et al.lll985t IGorenstein et al.lll988t ISahal 
2000; . Including an external convergence k would rescale 
the lens models shears by (1 — k). We consider how this 
would affect our comparison of model and environment 
shears in § 16.61 

As a fiducial case, we use all the available constraints, 
and we allow the lens galaxy's power law index to be a 
free parameter with a uniform prior in the range 0.2 < 
a < 1.8. Table [3] gives the range of recovered parameters, 
including the external shear, along with the \ 2 goodness 
of fit statistics for the best models. The \ 2 values are not 
always comparable to the number of degrees of freedom, 
particularly for WFI2033, for which x 2 = 52.06 for 5 de- 
grees of freedom. This suggests that there are subtleties 
in real lens galaxies that are not captured in our models 
(such as flux perturbations from substructure). Never- 
theless, we believe that the inferred shear distributions 
are reliable because omitting flux ratio constraints causes 
the shear distributions to shift by less than 2-sigma in 
all systems except RXJ1131 (which is affec ted by mi- 
crolensing and perhaps millilensing as well; iSluse et al. 
[2006|lChartas et al.ll2009l ; lCongdon et al.ll2010HDai et al. 
120101) . We have also tried imposing a prior on the lens 
galaxy's power law index to favor isothermal profiles. 

7 In Bayesian language, we adopt uniform priors; note that we 
work with quasi-Cartesian coordinates for the two components of 
shear, defined in Equation [T] below, and analogous coordinates for 
the two components of ellipticity. 
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None of these systematic effects (including flux ratios in 
RXJ1131) changes our conclusions about whether lens 
model shears do or do not match the shears estimated 
from the observed environments, so for simplicity we re- 
port only the fiducial results. The average shears in these 
six systems range from 7 = 0.06 to 0.28 for these models. 

4. ENVIRONMENT ANALYSIS: AN INDEPENDENT 
SHEAR DETERMINATION 

As a counterpoint to the shear inferred from lens mod- 
els, we want to directly calculate the shear due to exter- 
nal perturbations, including both the local group envi- 
ronment and the galaxies along the line of sight to the 
lens. Our main analysis tool is Monte Carlo simulations. 
We generate many possible realizations of the local lens 
environment and LOS such that the shears from the dif- 
ferent trials span the full range of shears allowed by our 
data. 

4.1. Overview 

Our goal is to use our spectroscopic and photometric 
data to constrain the mass distribution along the line of 
sight to each lens system in order to directly determine 
the external shear at the position of the lens. We de- 
scribe our formalism for calculating the external shear in 
§ 14.21 Several issues complicate this analysis. First of 
all, our spectroscopy is not complete down to our spec- 
troscopic magnitude limit, so we need to assign redshifts 
to those objects for which we do not have spectra. We 
assign redshifts based on the redshift distribution of the 
spectroscopically observed galaxies of similar I magni- 
tudes in each field. These galaxies that lack redshifts 
may include objects that are members of the lens host 
group, so we need to correct for incompleteness in the 
group membership as well. We detail our methodology 
for handling these incompletenesses in § 14.31 In addition, 
there are a number of observational and theoretical un- 
certainties (§ 14.41) that we consider in our analysis. To 
span a large range in parameter space including the spec- 
troscopic incompleteness and these uncertainties, we run 
Monte Carlo simulations to generate multiple realizations 
of the lens environments. 

There are two distinct components of the lens envi- 
ronment that we consider: the local group environment 
(§ I5.ip and the LOS galaxies that are projected close 
to the lens but that lie at different redshifts (§ I5.2[) . 
We assume that the group halo h as an NFW profile 
(|Navarro et al.l 119961 see Appendix IA.1I) and that the 
individual galaxies are trun cated singular isothermal 
spheres (TSIS; see Appendix IA.2[) . The galaxy velocity 
dispersions used in our shear calculati on are determined 
via t he Faber- Jackson relation (FJ; iFaber k, Jacksonl 
119761 ). We apportion mass between the group halo and 
the group galaxies to account for the varying degrees 
to which group members could be tidally stripped due 
to galaxy-galaxy encounters (§ I5.1.3[) . We apply cuts on 
projected separation from the lens and on apparent mag- 
nitude for the LOS galaxies (see ? 15.21 Appendix [Ct . We 
then calculate the shear from the full lens environment, 
including both the local group environment and the LOS 
galaxies (§ I5.3|) . 

For each Monte Carlo trial, we correct for spectroscopic 
incompleteness, apply our observational and theoretical 



uncertainties, and calculate the shear at the position of 
the lens for that particular realization. We run 1000 tri- 
als for each lens, and the resulting shear distributions 
represent the possible environmental shears allowed by 
our data. These distributions are then compared to the 
shears derived from lens modeling (§ [3]) to determine 
whether or not there is agreement between the two inde- 
pendent methods. 

4.2. Shear Determination 

The full lens potential contains contributions not only 
from the main lens galaxy, but also from all structures 
along the line of sight. For any perturber whose pro- 
jected offset from the lens is larger than the Einstein ra- 
dius of the lens, Re, we can expand the lens potential in a 
Taylor series and quantify the effect of the perturbation 
by the lowest-order significant terms in the expansion, 
convergence (k) and shear (7). k is a scalar quantity that 
represents the surface mass density of the system, but 
it cannot be inferred from image posi tions and flux ra- 
tios due to the mass-sheet degener acy (jFalco et al.lll985t 
iGorenstein et al.lll988HSahall2000D . The shear represents 
tidal distortions in the lensed image and can be expressed 
in terms of an amplitude 7 and direction # 7 (measured 
North through East). For a spherical perturber, 7 rep- 
resents the position angle of the perturber relative to the 
lens. We can also express the shear in terms of the two 
components 

7 c = 7cos(20 7 ) and 7 S = 7 sin(20 7 ). (1) 

The total shear, 7, is their quadrature sum, 7 = 
VtI+tF- It is useful to express the shear in these terms 
because shears from multiple perturbers at a given red- 
shift add linearly in (j c , 7 S ) space, whereas the shear am- 
plitude 7 does not. 

To determine the shear due to the full lens environ- 
ment, we use the multi-plane lens equation (Schneider 
et al. 1992; Petters et al. 2001) to handle perturbers 
at different redshifts. We discuss the full multi-plane 
shear formalism in a forthcoming paper (Keeton et al. in 
preparation) and highlight the key concepts here. The 
total shear is not simply the sum of contributions from 
individual lens planes; there are non-linearities because 
each perturber acts on light rays that have already been 
distorted by perturbers at other redshifts. To see this, 
consider the "shear tensor" for a perturber i, 



d 2 4>i 

().!■ 2 



dxidyi 
dxidyi dy'f 



Ki + 7c 
I si 



Jsi 
Hi - 1c 



(2) 



The shear tensors can be combined into the Jacobian 
matrix for the mapping between coordinates on the sky 
and coordinates in the source plane through the recursion 
relation 

j-i 

A$=I-£>ir 4 A<. (3) 

i=l 

where = I (the 2x2 identity matrix), and the sum 
runs over all perturbers (but does not include the main 
lens galaxy). Also, the weight factor /3y is a dimension- 
less combination of angular diameter distances between 
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TABLE 3 
Lens Model Results 



Lens 


tiE \ ) 


6 


PA (°\ 
FA ( ) 


Q 


7c 




A ,2 / nop 


HE0435 

PG1115 

RXJ1131 

HST14113 

B1422 

WFI2033 


1 9nl +u.uua 

±.zui_ 002 

1 1QC+0.004 
±.ioo_ 003 

l-889±g;S| 
0.839lg;gU 

n 7fi K+0.023 

u. ioo_ Q 0lg 
1 109+ 004 


u - iM -0.07 

0-261^4 
0.081°,°? 

u - Ji -0.06 

36+ 010 
U - 3D -0.12 


-ll.7t'( A l 

-83.7i*;° 
-«.o±?;« 

-56.0ti'g 
-87.3^ 


u -'^-0.24 
n c-,+0.16 
U-01_ .i2 

1 cq+0.14 

1 ic+0.23 
J — lo -0.28 

Q4+015 
u - 9 ^-0.13 
74+0-34 
u -'^-0.19 


0.064+^9 

-0.068t°;™ 

-0 n57+ 013 
u.uo(_ 020 

n 979 +0.067 
u - zlz -0.086 

-0 047+ - 013 
u - u *' -0.013 

247+ 018 
-0.022 


-0.038+roiS 

0.149±° : °i? 

001+ 008 
u.uu±_ 012 

o niQ+ 027 
u.uiy_ 041 

n i 79 +0.029 
u - 1,z -0.029 

1S6+ 0010 
u.iou_ 024 


20.49 / 5 

28.70 / 5 

27.71 / 5 
1.22 / 5 
2.22 / 5 

52.06 / 5 



Note. — For each quantity we quote the median and 68% confidence interval. The position angle is measured North through East. The 
power law index is defined such that the mass projected within radius R scales as M(R) oc R a 



the observer, planes i and j, and the source: 

DuD, 



D j D i 



(4) 



Note that /3y = 1 when j corresponds to the source 
plane. The total shear tensor is then defined by 



I — A' 



(5) 



If all shears are small (i.e., the components of Ti are 
< 1), we can make a Taylor series expansion and work 
to first order to approximate the total shear as 



(6) 



We account for all of the non-linear effects by con- 
structing the full multi-plane lens equation with each 
perturber along the line of sight included as a complete 
mass component (calibrated via the methods discussed 
in § EJ ■ This approach ensures that we properly handle 
both large shears from objects projected near the LOS 
and cumulative effects from galaxies at larger projected 
distances. We use the complete multi-plane formalism to 
compute the shear from the full lens environment. We 
also compare that with the shear obtained from the local 
lens environment alone to quantify the LOS effects on 
the shear. 

The definition of the shear we have adopted for our 
environment analysis is chosen because it enters the lens 
equation in the same way as the shear inferred from lens 
models, and we want the two quantities to be compara- 
ble. In principle, the lens equation can be manipulated to 
obtain other "flavors" of shear, which differ in the treat- 
ment of perturbers that lie outside the lens plane (details 
are given in Keeton et al. (in preparation) Although 
there are formal differences among flavors of shear, we 
find that they lead only to small shifts in the environment 
shear distribution and do not change the conclusions we 
draw when comparing environment shears to lens model 
shears. The insensitivity of our results to the shear fla- 
vor may arise because we have selected lenses that lie 
in group environments, which have a large component of 
shear from the main lens plane. 

All of our shear calculations are performed using an 
updated version of the software developed by iKeetonl 



8 For example, Momcheva ct al. (2006) use the "effective shear", 
which gives reduced weight to perturbers at higher redshift offsets 
from the main lens. 



(|2001l) . which takes as inputs the relative positions of 
the lens and perturbing masses, as well as their Einstein 
radii, redshifts, and truncation radii. The code returns 
7 C and 7 S for the full input environment and for the local 
group environment alone, the latter including only ob- 
jects within a small redshift Sz of the lens plane. We set 
the redshifts of all lens plane galaxies, including the lens 
itself, to the mean redshift of all group member galaxies 
to eliminate redshift-space distortions due to the peculiar 
velocities of the group galaxies. 

4.3. Correcting for Spectroscopic Incompleteness 

We need to know the redshifts of the galaxies in our 
sample to properly compute the environment contribu- 
tion to the shear. Our spectroscopic data do not fully 
sample all of the objects in our photometric catalogs 
down to the spectroscopic limit of I = 21.5 (Figure [3]), 
affecting our measurement of both the local lens envi- 
ronment and the line of sight. There are likely group 
members without spectra (i.e., those that lie on the red 
sequence at the lens redshift), so we correct for these in 
considering the local environment. We also assign red- 
shifts to the likely LOS perturbers without spectroscopy 
based on the redshift distribution of our spectroscopic 
sample as a function of apparent total magnitude. For 
each Monte Carlo trial, the local group environment in- 
cludes the group dark matter halo, the spectroscopically 
confirmed group members, and the photometric red se- 
quence galaxies that are assigned group membership in 
that trial. The LOS perturbers include the spectroscopi- 
cally confirmed non-group members in the field, the pho- 
tometric red sequence galaxies not assigned group mem- 
bership in that trial, and the non-red sequence galaxies 
without redshifts. 

To correct for spectroscopic incompleteness in the 
host group, we consider galaxies without redshifts that 
lie on the photometric red sequence at the lens red- 
shift and assume that some are group members. The 
"red sequence" galax ies used here and described by 
iWilliams et al.l (|2006t ) and Williams et al. (in prepa- 
ration) are defined within a very narrow range of col- 
ors, roughly within ±0.05 magnitudes, whic h is the typ- 
ical l cr scatter of the red sequence (e.g. iBower et al.l 
Il992t iMcIntosh et al.l 120051 ). This definition is conser- 
vative in that it works to exclude interlopers at other 
redshifts, but it also excludes some true group members 
with old stellar populations that do not quite fall within 
this narrow color band. All galaxies, red or blue, that 
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are spectroscopically-confirmed group members are still 
accounted for in our analysis of the lens environment. 
We assign each of those galaxies a probability of group 
membership equal to the fraction of spectroscopically ob- 
served red sequence galaxies that are confirmed group 
members (see Table [2]). For each Monte Carlo realiza- 
tion of the environment, every red sequence galaxy has 
this probability of being chosen for group membership 
in that trial. If a particular red sequence galaxy is not 
assigned group membership in a trial, we classify it as 
a LOS galaxy. For different Monte Carlo trials, the red 
sequence galaxies assigned group membership can there- 
fore vary in number, although every individual galaxy 
has the same probability of being added to the group as 
any other red sequence galaxies in a trial. In three cases 
(HE0435, Q0047, WFI2033), no red sequence is found 
near the lens redshift, despite the presence of a spectro- 
scopically confirmed group. For these groups, we do not 
assign membership to any objects that are not already 
confirmed group members. 

We assign redshifts to LOS galaxies without spectro- 
scopic data based on the redshift distribution of our 
spectroscopic sample in that particular field, excluding 
galaxies confirmed to be in the host group of the lens. 
Rather than combine the redshift distributions for all 
systems into an aggregate catalog, we create separate 
distributions to preserve the large-scale structure in red- 
shift space in each field. We restrict our spectroscopic 
sample to objects projected within 5' of the lens and 
brighter than I = 21.5 (little shear is contributed by ob- 
jects outside these limits, see Appendix ICl) . We separate 
the galaxy redshift distributions into 8 bins of total I- 
band magnitude on a system-by-system basis. Each bin 
has a width of ^0.5 mag except for the first bin, which 
includes all objects brighter than / ~ 18. The width and 
locations of the bins are adjusted slightly to include at 
least three objects so that we can discriminate between 
over- and underdense regions along the line of sight. The 
bin centers and widths rarely vary by more than 0.2 mag- 
nitudes in either direction. For each bin in each system, 
we sort the objects by redshift and compute the cumula- 
tive probability distribution for redshift. For each object 
in the photometric catalogs without spectroscopic data, 
we draw a random number distributed uniformly in the 
range [0,1] and assign it the redshift at which the cumu- 
lative distribution for the corresponding magnitude bin 
attains that value. 

We make a slight modification when applying this red- 
shift assignment procedure to photometric red sequence 
galaxies that are not selected for group membership. For 
these galaxies, we again draw from a redshift distribution 
in bins of total /-band magnitude as described above. 
However, because red sequence interlopers tend to be in- 
trinsically blue galaxies at higher redshifts than the lens, 
we use parent distributions including only spectroscopic 
galaxies with z > zi ens . This limits the number of objects 
in each magnitude bin, particularly for the bright (/ ~ 
18-19) galaxies. To maintain our condition of having at 
least three objects in each bin, we increase the faint-end 
limit of the brightest bin to / ~ 19 and restrict ourselves 
to 6 bins instead of 8. 

One caveat with this method is that the probability of 
selecting a redshift between two adjacent perturbers in 
redshift space is assumed to be uniform over that range. 



This approximation is fine for well-sampled overdense re- 
gions, but does not handle voids in redshift space prop- 
erly. Overdense regions generally induce larger pertur- 
bations on the lens potential than underdense regions, so 
we are less concerned with voids (see § I6.5[) . 

4.4. Additional Sources of Error 

Once all redshifts are specified, we use our observa- 
tional data to build a model of the mass distribution 
along the line of sight. There are uncertainties in the ob- 
servational data themselves that could affect the shear, 
including: 

• Slope and scatter in the Faber- Jackson relation, 
which we use to convert luminosities into veloc- 
ity dispersions and Einstein radii. We allow for an 
error of 0.20 in the slope, as well as an intrinsi c 
scatter of ~0.07 in loga (jBernardi et al.l l2003al ) . 
These errors affect all galaxies in both the local 
group environment and along the line of sight. 

• Error in the group centroid position, which sets 
the position of the group dark matter halo. In 
our shear analysis, we generate a centroid position 
from a bootstrap resampling of the group galaxy 
positions with uniform weights. The group cen- 
troid errors are determined from the standard de- 
viation of 1000 bootstrap trials (Table [2]). Using 
luminosity- weighted centroids does not change our 
conclusions. 

• Error in the group velocity dispersion, which de- 
termines the total group mass. We use a boot- 
strap resampling of the group g alaxy redshifts and 
a bi- weight estimator for scale (|Beers et al.lll990l) 
to determine the group velocity dispersion and er- 
ror (Tabled]). 

• Photometric errors in the galaxy magnitudes 
(Williams et al. in preparation), which we use to 
calibrate the mass models for individual galaxies in 
both the local group environment and the LOS. 

There are also uncertainties associated with the theo- 
retical assumptions we must make: 

• Choice of density profile for the group dark matter 
halo. We adopt NFW profiles for our fiducial mod- 
els, although we also consider a singular isother- 
mal sphere (SIS). Details of the difference between 
SIS an d NF W group halos are presented in Ap- 
pendix IA.1I 

• Scatter in the concentration parameter, c V i r , of 
the group dark matter halo. We allow for a scat- 
ter of 0.14 in \ogc V i r , which is approximately con- 
stant over a mass range that encompasse s the virial 
masses of the groups in ou r sample (Bullock et al.l 
[2001 IWechsler et al.l[200l . 

• The apportionment of mass between the group halo 
and the individual group galaxies (§ 15.1.31) . 

We build all of these uncertainties into our Monte Carlo 
analysis by drawing values from the appropriate distribu- 
tions for each trial. We examine the relative importance 
of the various uncertainties in our results in § 16.11 
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5. SHEAR FROM ENVIRONMENT COMPONENTS 

5.1. Shear Contribution from Group Environment 

The local group environment of each lens consists of a 
group dark matter halo and the individual group galax- 
ies. In principle, the group can contain substructure 
on a smaller scale than the galaxies brighter than our 
spectroscopic limit, but for this analysis, we will assume 
that its contribution is small in comparison to the overall 
group halo. We only include spectroscopically-confirmed 
or Monte Carlo-assigned group members within a virial 
radius of the group centroid in each Monte Carlo trial 
because the virial mass of the group should not be ap- 
portioned to objects outside this radius. We treat group 
galaxies outside the virial radius as LOS galaxies at the 
lens redshift. 

5.1.1. Shear from Group Halo 

In modeling the group dark matter halo, it is not 
clear whether virialized groups have mass density pro- 
files that are better fit b y a singular isother mal sphere 
(SIS) or a NFW profile ((Navarro et al.1 119961 ). Further- 
more, if the group is not yet virialized, the form of the 
halo profile could be very different from either profile. 
This possibility is somewhat mitigated by the fact that 
non-virialized groups are likely to have a larger fraction 
of their total mass in individual group galaxies, which 
may not yet have b een stripped via tidal intera ctions 
with other galaxies (jZabludoff fc Mulchaevl [19981. We 
test how much of an effect the choice of an SIS or NFW 
profile has on the shear from the group halo and find 
that in the majority of cases, the effect is small, particu- 
larly when much of the group mass is apportioned to the 
individual galaxies rather than the overall dark matter 
halo. The d etails of the shear calculations are presented 
in Appendix lA.il We assume an NFW profile throughout 
our analysis because the choice of group profile does not 
significantly affect on our results, and NF W profiles are 
bette r fits to group-size dark matter halos (|Navarro et al.l 
fl996l) . 

5.1.2. Shear from Individual Group Galaxies 

We assume that each group galaxy is a SIS in the cen- 
tral regions, which has been shown to be an accurate 
model for galaxy halos (e.g. iRusin fc Kochanekl 120051: 
Koopmans et al.ll2006t iGavazzi et al.ll2007b iBolton et al.l 
20081: iNipoti et al.ll2008t iKoopmans et al.ll2009D. Realis- 
tic g alaxy are truncated (see also e.g. ISuvu fe Halkolal 
2010), so we assume that the profiles are truncated 
at radius r t , as described in § 15.1.31 and Appendix [Bj 
The shear calculation for a truncated singular isothermal 
sphere (TSIS) depends on the internal velocity dispersion 
cr, from which we can determine the Einstein radius, Re- 
The de tails of this shear calculation are presented in Ap- 
pendix E2J 

For each trial, the velocity dispersions of the confirmed 
or assigned group galaxies are determined from their ob- 
served total /-band magnitudes via the Faber - Jackson re- 
lation (hereafter FJ; iFaber fc Jackson! 1 19761 ). L cx a lFJ . 
This procedure assumes that all group members lie on 
this relation (i.e. have kinematics like those of early-type 
galaxies). As a result, we are effectively assigning up- 
per limits to the masses of any late type galaxies, which 
tend to have lower masses than what would be calculated 



from this relation. This assumption is crude, but neces- 
sary, given that we currently lack detailed morphologies 
for the galaxy sample. It is not baseless, as we estimate 
that roughly ~ 60 — 80% of the galaxies in the host group 

are red and/or non-star forming (see § I6.5fl . 

We a dopt the FJ relation given bv iBernardi et ahl 
(12003a!) for the SDSS r-band, 



logo- = 2.2 



0A(M r + 21.15) 



3.91 



^slope 



(7) 



where e s i ope and Cj nt are the FJ sl o pe erro r and intrinsic 
scatter term from [Bcrnardi et al.l (|2003aD , respectively. 
These terms are assumed to be Gaussian distributed with 
a rms equal to 0.20 and approximately 0.07 respectively 
for the r-band. Of these two terms, the intrinsic scat- 
ter term dominates the effect on our resulting shear dis- 
tributions, so future references to the effect of the FJ 
uncertainty on the shear reflect primarily the scatter. 

We calculate the K-correction between the Cousins I- 
band and the SDSS r-band by synthesizing colors from 
an elliptical galaxy template derived using stella r popu - 
lation synthesis models from iBruzual fc Chariot! {2003) . 
We assume a passively evolving simple stellar popula- 
tion with solar metallicity and a Salpeter IMF formed at 
z = 3 and evolved to the present epoch. If the majority 
of our group galaxies are ellipticals, this shou ld provide a 
valid approximation. IBernardi et al.1 (|2003a| ) also include 
a passive evolution term in the determination of M r that 
is linear in redshift. We omit this te rm and account for 
passiv e evolution explicitly using the IBruzual fc Chariot! 
(2003) elliptical galaxy model because the range of red- 
shifts over which we are correcting is larger than for the 
SDSS sample, where a linear correction is a better ap- 
proximation. For each galaxy, we evolve the model from 
the galaxy's redshift to z = 0.1, which is roughly the 
median redshift of the SDSS sample. Our determination 
of M r for a particular galaxy at redshift z ga i is 

M r = I obs - DM{z ga i) +K Itr (z gal ) + E I (z gal , 0.1), (8) 

where DM{z ga i) is the distance modulus to the galaxy, 
K^ r is the K-correction term, and Ei is the passive 
evolution term . We are unable to use the Fundamen- 
tal Plane (e.g. [Diorgovski fc pavisj|1987l; IBernardi et al.l 
!2003bl : !Cappellari et al.ll2006HRobertson et al.ll2006D be 
cause we currently do not have effective radii for the 
galaxies. 

From this relation, we calculate the Einstein radius, 
Re, for each group galaxy. Re for a general SIS per- 
turber along the line of sight is 



Re = 



Air a 2 Dps 
c 2 D S ' 



(9) 



Re depends on the angular diameter distance to the 
source (Ds) and the angular diameter distance from the 
perturber to the source {Dps)- For the group galaxies, 
Dps is equal to D^s, the angular diameter distance from 
the lens redshift to the source. 

The Einstein radius of each lens galaxy is approxi- 
mated as half the imag e separation, obtain ed from the 
CASTLeS database and IRusin eraTI (|2003| ). These val- 
ues are within a few hundredths of an arcsecond of those 
determined from lens models of the 4-image lenses in 
our sample, so this is a good approximation. We do not 
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use the FJ relation for the lens galaxies themselves be- 
cause our photometry is not always able to accurately 
deblend the galaxy from the lensed images. We also do 
not use the lens galaxy fluxes from the CASTLeS web- 
page due to an unexplained discrepancy between those 
magnitudes and ours in the few systems where we can 
deblend our images. In those systems, the velocity dis- 
persions determined from applying the FJ relation to our 
magnitudes are consistent with both those derived from 
assuming Re as half the image separation and from di- 
rect spectroscopic measurements of the internal galaxy 
kinematics. This is not true for the CASTLeS magni- 
tudes, which generally are anomalously faint. The prop- 
erties of the lens galaxy itself do not directly enter into 
the calculation of the external shear, and are only used 
in determining the mass distribution within the group 
(§ I5.1.3|) . For these reasons, we treat the lens galaxies in 
this separate manner. 

5.1.3. Reapportioning Mass between Group Halo and 
Individual Galaxies 

iMomcheva et al.l (|2006l ) calculate the expected shear 
for lenses in six galaxy groups in the limits where the 
shear is either due to a group halo or to galaxies at the 
group member positions. In our analysis, we examine 
the more realistic case where there is mass in both the 
group halo and the galaxies by redistributing the total 
group mass between the two components. 

We define the parameter fhaio to be the ratio of the 
group ha lo m ass to the virial mass of the group (see 
Equation EU, 

halo = (10) 

We expect galaxy halos to be tidally stripped as the 
group forms and ages. In the extreme limit fhaio = 0, 
all of the mass is attached to the galaxies and there is 
no mass in a common halo. As we approach large fhaio, 
virtually all the mass has been stripped from the galaxies 
and merged into a common halo. Therefore, fhaio char- 
acterizes the extent to which the halos of group galaxies 
may have been tidally stripped. Because we do not have 
observational constraints on fhaio, we allow it to vary as 
a free parameter in our Monte Carlo simulations. 

For each realization of the group environment in our 
Monte Carlo simulations, we choose a random value of 
fhaio uniformly distributed between and 1 and assign 
masses to the group halo and the group galaxies as de- 
scribed above. We distribute the total mass, M vir , be- 
tween the group halo and the confirmed or assigned group 
members projected within r„j r of the group centroid. For 
a particular value of fhaio, the galaxies' truncation radii, 
Tt, are scaled so that the density at the truncation radius 
is the same for all group galaxies, and the total mass of 
the galaxies is equal to (1 — fhaio)Mvir- The details are 
in Appendix |B| 

Since the mass apportionment can vary depending on 
the distribution of galaxy luminosities for the group 
members, the values of fhaio are not directly compara- 
ble from group to group in terms of how fhaio affects 
the shear. Similar values of fhaio between groups do not 
necessarily imply that an L* galaxy has the same trunca- 
tion radius from group to group. For the group with the 
most mass per galaxy in our sample (HE0435), a value of 



fhaio — implies typical truncation radii of over ~1 Mpc, 
while for the least massive group (MG1654), it implies 
truncation radii of ^60 kpc. 

For large values of fhaio, the truncation radii of the 
group galaxies may be smaller than their typical lumi- 
nous radii, which is unphysical. For Monte Carlo tri- 
als where this occurs, we repeat the trial with a new 
value of fhaio unti l the trunca t ion ra dii become phys- 
ical. A study by Ro ssi et all (|2010l ) using early-type 
galaxies from the SDSS finds that the effective (half- 
light) radii of galaxies with absolute magnitudes simi- 
lar to those of our galaxies are best fit by the relation 
log (.R e //) = -0.257M r - 5.086. We determine each 
group galaxy's absolute r-band magnitude, use this re- 
lation to calculate R e ff, and take 2i? e // to be the min- 
imum truncation radius allowed. For a de Vaucouleurs 
profile, this radius encloses ~ 70% of the total light of 
the galaxy. A Monte Carlo trial is repeated with a new, 
smaller value of fhaio if any galaxy in the trial is assigned 
a truncation radius smaller than this. We run tests in 
which the minimum allowed truncation radius is 5i? e // 
(enclosing ~ 90% of the light), but this does not change 
our results. 

5.2. Shear Contribution from Line-of-Sight Objects 

Objects along the line of sight to the lens also con- 
tribute to the extern al shear by in ducing perturbations 
in the lens potential (jKeetonll2003l ). even if they are not 
dynamically associated with the group. Our shear cal- 
culation software can determine the contribution from a 
perturber at any redshift along the line of sight, making 
it straightforward to calculate the contributions from all 
LOS objects. We apply a similar formalism to the LOS 
galaxies as we do for the group galaxies (§ 15.1. 2j) . Using 
Equations [7] and we determine their velocity disper- 
sions (and therefore, their Einstein radii) from their ob- 
served total magnitudes and redshifts. Our use of the FJ 
relation is potentially more problematic here, but we still 
estimate that ~ 50% of the LOS galaxies are red and/or 
non-star forming (see § 16. 5|) . The truncation radii of the 
assumed TSIS for the LOS galaxies are approximated as 
^200 ! where the mean density inside the volume bounded 
by a sphere of that radius is 200 times the mean matter- 
density of the universe. r2oo for an isothermal sphere at 
a redshift z ga i and with a velocity dispersion a is 



10H oy /Q m (l + z ga i) 3 ' 

To account for all LOS perturbers, we need to correct 
for the spectroscopic incompleteness for objects brighter 
than our limit of I = 21.5. Our photometric catalog con- 
tains many galaxies brighter than this limit without spec- 
troscopic data. We model their redshift distribution in a 
physically motivated way (§ 14. 3p . calculate their masses 
and truncation radii (Equation II 1 [) , and determine their 
effects on the lens shear. 

One concern is the effect of LOS objects projected far 
enough away from the lens that they lie outside of our 
sampling region. Including objects projected far from the 
lens is both observationa lly and computa tionally unfea- 
sible. Simulations run by iBrainerdl (|2010f) show that the 
external shear due to galaxy-galaxy lensing extrapolates 
to zero at a projected separation of ~5'. Observational 
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studies (e.g. lAuger et al.ll2007l ) have suggested that the 
most significant perturbers are projected even closer to 
the lens (~15"). Given that the shear contribution from 
an object decreases with increasing project ed sep aration, 
we make and justify a cut at 5' fAppendix lC.l[) . 

Finally, we consider the shear contribution from ob- 
jects fainter than our spectroscopic limit, which is more 
difficult because we do not have a large spectroscopic 
sample at fainter magnitudes from which to build redshift 
distributions. To avoid this problem, we cut our photo- 
metric sample at our spectroscopic limit oil = 21.5. We 
gauge the effects of faint er ob jects and give justification 
for this cut in Appendix IC.21 We find that a cut at / = 
21.5 includes the most significant perturbations and is a 
reasonable cutoff magnitude for our analysis. 

5.3. Shear from Full Lens Environment 

The local group environment and LOS perturbers both 
contribute to the observed shear, although by varying 
amounts. Whereas the previous sections outlined our 
methodology for determining their respective shear con- 
tributions, a comprehensive treatment of the shear cal- 
culation must simultaneously include all of the pertur- 
bations from the full lens environment. Considering the 
LOS and local environment effects separately may be ap- 
proximately correct if the effect of one component is small 
relative to the other. However, we do not know a priori 
whether this is the case, and if so, which component is 
the dominant one. Such an analysis does not account for 
interactions among multiple lensing planes, so we run a 
full simulation with all perturbers included in order to 
properly handle these effects. At a fundamental level, 
higher-order terms in the lens potential can be analyzed 
in the full analysis (Keeton et al. in preparation), but 
we only consider the shear term in our calculation. 

In our Monte Carlo simulations of the full lens environ- 
ments, we run 1000 trials for each lens. For each trial, 
we assign group membership to a subset of photomet- 
ric red sequence galaxies based on the fraction observed 
from spectroscopy (S I4.3[) and draw randomly from a dis- 
tribution of fhaio values with a uniform prior (requir- 
ing that no galaxies are assigned a truncation radius less 
than 2R e ff). The total group mass is then apportioned 
among the group halo and the individual group members 
(§ 1 5 . 1 3 [) . We assign random redshifts to the LOS galax- 
ies using the procedure in § 14.31 We assume Gaussian 
errors to account for most of the uncertainties discussed 
in § 14.41 except when stated otherwise. Each trial there- 
fore represents one possible realization of the full lens 
environment. By running 1000 trials, we are sampling a 
large region of our parameter space in order to determine 
the range of external shears that could be produced by 
the environment. 

6. RESULTS & DISCUSSION 

6.1. Effects of Errors on Full Lens Environment Shear 

We run five sets of 1000 realizations of the full lens 
environments (§ I5.3j) . allowing for variations due to one 
of our uncertainties (§ I4.4p in each set. We also run a set 
of trials representing a control sample excluding these 
uncertainties, as well as a set of trials including all of 
the uncertainties simultaneously (Figure [4j. The varia- 
tion in the shears for the control set essentially accounts 



just for variations in fhaio, the selection of additional 
group members from the photometric red sequence, and 
the variations in the redshift distributions of the LOS 
perturbers, but for none of the additional uncertainties 
described in § 14.41 

The mean and scatter in 7 C and 7 S for each uncertainty 
is listed in Table |U We assume Gaussian errors in cal- 
culating the errors in 7 C and 7 S , but the actual shear 
distributions can be asymmetric and non-Gaussian (see 
Figured]). We note that: 

• The scatter in the shear for the control sample is 
typically ~0.01, whereas the scatter when all uncer- 
tainties are taken into account is typically ^0.03. 
The average increase in the scatter in either 7 C or j s 
is ^0.015, although this varies among the individ- 
ual systems. The means of the shear distributions 
themselves have a typical shift of ^0.01, although 
7 C in HE0435 changes by 0.05, the largest offset for 
any system. 

• The scatter in the Faber- Jackson relation and un- 
certainty in the group centroid position are usually 
the dominant sources of error. The effect of the FJ 
scatter on the shear is generally more significant for 
low values of fhaio, when more mass is assigned to 
the group galaxies and their shear contribution is 
more greatly affected. Conversely, the uncertainty 
in the group centroid position is generally more im- 
portant for high values of fhaio, as the group halo 
has a larger effect on the shear when it is more 
massive. 

• The magnitude errors in the photometry are neg- 
ligible, resulting in changes in the mean and stan- 
dard deviation of 7 C or 7 S of less than 0.01. 

• The offset and increase in scatter in the shear due 
to the group velocity dispersion error and concen- 
tration parameter, both of which affect the group 
mass, are on the order of < 0.01. 

• There is a large amount of scatter with little sys- 
tematic dependence on fhaio- As a result, we can- 
not place useful constraints on fhaio by comparing 
shears from lens model results to the shears from 
the environment (§ I6.2[) at this time. 

• For the full lens environment, the average shear is 
(7) = 0.08, ranging from as little as 0.02 to as much 
as 0.17. The environmental effects on the lens po- 
tential are therefore comparable to the lens model 
shears (Table [3]) and cannot be ignored when an- 
alyzing the lens potential. We consider the shears 
due to the local group environment alone in § 16.31 

We also run trials where we take into account only 
the galaxies for which we have spectroscopic data. Our 
results are qualitatively similar except in the cases of 
B1422 and WFI2033, two of the systems where our shears 
do not match the lens model results (§ 16. 2[) . For these 
two systems, the shear distributions shift substantially 
but in a direction that makes them more discrepant with 
the lens model shears. 
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Fig. 4. — 7 C vs. 7„ for the full lens environment with only systematic effects (variations in fhalo? selection of additional group member 
candidates from the red sequence, variations in the redshift distribution of the line-of-sight perturbers; top row), isolated observational 
uncertainties (2nd through 6th rows), and the systematics plus all observational uncertainties (bottom row). The points arc color coded 
by fhalo as indicated by the color bar. For the six 4-image lenses, the la and 2a error ellipses derived from lens modeling are shown for 
power-law lens galaxy models (bottom row). The scatter in each panel is quantified in Table|4] In three of the six 4-image lenses (HE0435, 
PG1115, HST14113), our environmental shears match t he m odel shears, but in the others (RXJ1131, B1422, WFI2033), the environmental 
and model shears are inconsistent at the > 95% level (§ 16.21 1. 
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TABLE 4 

Shear Statistics for Observational Uncertainties 
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Note. — Errors are standard deviations of the shear distributions assuming a Gaussian distribution. In some cases, the Gaussian 
assumption is very rough (see Figure 0J, and these errors should be considered as guides only. 



6.2. Comparison of Environment and Lens Model 
Shears 

Our shear distributions overlap the lens model-derived 
values in HE0435, PG1115, and HST14113, but not in 
RXJ1131, B1422, or WFI2033 (Figure g]). We want to 
quantify the level of disagreement between the lens model 
and environment shears, which are determined indepen- 
dently from one another. To do this, we construct the 
probability distribution for A7, the offset between the 
shear distributions, to determine whether it is consistent 
with zero. If we have probability distributions p e nv{.l) 
and Pmod{l) for the environment and lens model shears 
respectively, the formal statement of the probability dis- 
tribution of A7 is 

p(A 7 ) = / Penv {l')Pmod{i + A 7 )d 7 '. (12) 

For the environment shear, we have discrete samples 
rather than a continuous distribution of shears. We can 
also sample 1000 discrete points in ( 7c , 7s ) space from 
the lens model shear distributions. We cross-correlate 
these points with the 1000 shear points from the en- 
vironment measurement and determine the difference 
between each model shear i and environment shear j, 
A 7ij = 7 ™ od - jf nv . The discrete estimator of the A 7 
probability distribution from these samples is then 

, N m r N mo d 

p(A 7 ) = „ „ E E <5(A 7 + 7 f"'- 7 r od )- (13) 

1=1 7 = 1 



In other words, we have a collection of N env N mo d = 10 6 
points at the locations A 7 y for all i and j in both 7c and 
7s separately. We plot these points in (A 7c ,A 7s ) space 
and calculate isoprobability contours. The P- value for 
A 7 = 0, which represents perfect agreement between the 
two distributions, is the fraction of points lying outside 
the contour passing through the origin. The 68%, 95%, 
and 99.7% contours are shown in Figure [5] and the value 
of -P(A 7 = 0) is given for each system. 

The results show that the environment and lens model 
shears for HE0435, PG1115, and HST14113 agree to 
within the 95% contours, albeit marginally for PG1115 
and HST14113. However, the shears for RXJ1131, 
B1422, and WFI2033 do not. The disagreements could 
arise from a number of effects that we have not consid- 
ered in our environment analysis, as well as problems in 
the lens modeling itself. These are discussed further in 
§E!and §EH 

6.3. Constraining Line-of-Sight Shear 

We compare the shear due to the local group envi- 
ronment alone to that of the full environment to quan- 
tify the effects of the LOS perturbers on the shear (Fig- 
ured]). The local lens environment alone induces an av- 
erage shear of 7 = 0.05 with a range from 0.01 to 0.14, 
compared to that for the full environment ( 7 = 0.08, 
ranging from 0.02 to 0.17; Table [S]). The shear ampli- 
tude from the LOS objects can thus be comparable to 
that from the local group environment, although it de- 
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Fig. 5. — 68%, 95%, and 99.7% contours for A7 in each of the six 4-image systems in our sample. The origin represents perfect agreement 
between the environment shears and the lens model shears. The P-value for A7 = is given for each system. A general power-law is 
assumed for the lens galaxy density profile in the lens modeling. For three of the systems (HE0435, PG1115, HST14113), the environment 
shears are consistent with the model shears to within the 95% confidence intervals. However, the shears are inconsistent at greater than 
95% confidence for the remaining three systems (RXJ1131, B1422, WFI2033). 



pends on the configuration of the individual lens systems. 
We also compute Aj env ^ the mean offset in {pf c , 7 S ) space 
between the shears from the local environment alone and 
the shears from the full environment. The mean value of 
this offset across all systems is A^ env = 0.06 with a range 
from 0.02 to 0.11, indicating that the LOS can be a signif- 
icant perturbation to the lens potential. Accounting for 
the LOS objects does not significantly add to the shear 
scatter by more than 0.02 in any system. 

6.4. Notes on Individual Lens Systems 

From the results in Figures |4] and [6] and quantified 
in Tables 2] and [5j it is clear that the behavior of the 
shear is not uniform across all lenses. This is expected, 
as group environments can vary significantly in overall 
mass, richness, and other properties that can affect the 
lensing potential. In addition, the configuration of indi- 
vidual galaxies relative to the lens galaxy can have large 
shear effects, particularly in the limit of low fhaio- The 
shear from LOS galaxies can be sensitive to a small num- 
ber of objects projected close to the lens, as well as to 
the large amount of cosmic varia nce in beams along lines 
of sight to individua l systems (|Momcheva et alj l2006t 
IFassnacht et al.l I2010f ) . In this subsection, we provide 
some qualitative insights about the effects of the group 
environments and LOS objects in individual systems, as 
well as previous lens model results. 

The scatter in the shear distribution for RXJ1131 is rel- 
atively small. One reason is that the lens galaxy is offset 
from the group centroid more than in any other group, so 



perturbations from the group halo and most of the other 
group galaxies have less of an impact on the lensing po- 
tential. There is evidence for subst ructure in RXJ1131 , 
resulting in anomalous tim e delays ([Morgan et aJ]l2006t 
iKeeton fc Moustakasl 120091 ICongdon et al.l 12010ft . The 
lens model shears become even more discrepant from our 
environment results if the flux ratios are not taken into 
account. 

Accounting for LOS objects in B1422 changes the shear 
by Aj env — 0.11, the largest shift for any system. Much 
of this is likely due to the presence of a bright (/ = 18.38), 
nearby LOS galaxy at a projected separation of 0.13'. 

The shear from the local environment of WFI2033 
alone is small, but accounting for LOS objects shifts the 
shear by Aj env — 0.08. This is likely due to the presence 
of a nearby LOS object projected 0.06' from the lens, as 
well as a few nearby LOS objects in our photometric cata- 
log, similar to B1422. There are six objects brighter than 
/ = 21.5 that are projected between 0.12' and 0.31' away 
from the lens that ar e not in our spectroscopic catalog. 
ICongdon et al.l (|2010f ) find time delay anomalies in this 
system, suggesting either a complex lens environment or 
the presence of substructure. 

The scatter in the shear distributions due to uncer- 
tainty in the position of the group centroid is the domi- 
nant source of uncertainty in HE0435 and PG1115, and 
also has a significant effect in Q0047 and HST14113. The 
panel in Figure [4] showing the effects of the centroid error 
indicate that for large values of fhaio, the shear ampli- 
tude is shifted a larger distance away from the typical 
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Fig. 6. — 7c and ■js from the lens group environment alone (left panels) and from the full environment including line-of-sight galaxies (right 
panels). The right panels are identical to the bottom row of Figure [4] Each point represents one of 1000 realizations of the environment, 
including the spectroscopic or photometrically assigned group galaxies, the LOS galaxies, and the uncertainties in Figure [4] The points 
are color coded by fhalo as indicated by the color bar. For the six 4-image lenses, the la and 2a error ellipses derived from lens modeling 
are shown for power-law lens galaxy models. Statistics for each distribution are in Table [5] Including the LOS galaxies has a large effect 
on the shear in B1422 and WFI2033 due to perturbations from LOS objects projected close to the lens. 
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TABLE 5 

Shear Statistics From Local Group Environment and Full 
Environment 



Lens 


i^c) Local 
(Ts) local 
local 


{lc)full 

(is) full 

<7> full 




Q0047 


0.02 ± 0.03 
0.00 ± 0.02 
0.03 ± 0.03 


0.04 ± 0.03 
0.01 ± 0.02 
0.05 ± 0.02 


0.05 ± 0.03 


HE0435 


0.02 ± 0.03 
-0.03 ± 0.04 
0.05 ± 0.03 


0.02 ± 0.04 
-0.01 ± 0.05 
0.05 ± 0.04 


0.07 ± 0.05 


MG0751 


-0.12 ± 0.05 
-0.05 ± 0.03 
0.14 ± 0.05 


-0.15 ± 0.05 
-0.08 ± 0.03 
0.17 ± 0.05 


0.08 ± 0.04 


PG1115 


-0.05 ± 0.04 
0.04 ± 0.03 
0.08 ± 0.03 


-0.07 ± 0.04 
0.06 ± 0.03 
0.10 ± 0.03 


0.06 ± 0.04 


RXJ1131 


-0.01 ± 0.01 
0.00 ± 0.005 
0.01 ± 0.01 


-0.01 ± 0.01 
-0.01 ± 0.01 
0.02 ± 0.01 


0.02 ± 0.01 


HST14113 


0.07 ± 0.04 
-0.03 ± 0.04 
0.09 ± 0.04 


0.11 ± 0.04 
0.02 ± 0.04 
0.12 ± 0.04 


0.09 ± 0.05 


B1422 


-0.02 ± 0.02 
-0.05 ± 0.03 
0.05 ± 0.03 


0.05 ± 0.02 
-0.13 ± 0.03 
0.14 ± 0.03 


0.11 ± 0.03 


MG1654 


-0.00 ± 0.01 
-0.00 ± 0.01 
0.01 ± 0.01 


-0.02 ± 0.01 
-0.01 ± 0.01 
0.03 ± 0.01 


0.03 ± 0.01 


WFI2033 


0.01 ± 0.01 
0.00 ± 0.01 
0.01 ± 0.01 


-0.03 ± 0.04 
0.06 ± 0.02 
0.08 ± 0.03 


0.08 ± 0.03 


Average 


0.05 


0.08 


0.06 



Note. — Errors are standard deviations of the shear distri- 
butions assuming a Gaussian distribution. In some cases, the 
Gaussian assumption is very rough (see Figure [Sj , and these er- 
rors should be considered as guides only. 



shears in the corresponding control samples. The sepa- 
ration between high and low values of fhaio is not clearly 
defined, so we cannot constrain fhalo based on a com- 
parison of our shear distribution to the lens model shear 
distribution. We might infer that there arc more real- 
izations of the environment with low fhaio that match 
the lens model-derived shears. Better spectroscopic sam- 
pling of the group members in this system would help 
to tighte n the centroid position errors and reduce the 
scatter. iKochanek et al.l ([2000; performed a modeling 
analysis of HE0435 assuming an ellipsoidal pseudo-Jaffe 
model for the lens, an SIS model for a nearby perturber 
(one of the objects in our photometric red sequence), plus 
external shear. They find a best-fit shear amplitude of 
7 ~ 0.05 and Q 1 ps —30°, which is consistent with our 
shears calculated from the environment when converted 
to (7c, 7 S ) space. However, their results are not directly 
comparable to our lens model results as we do not model 
individual pcrturbers. 

In the host group of HST14113, the brightest group 
galaxy (/ = 16.82) is located only 0.68' from the lens, 
indicating that it could have a large effect on the shear. 
However, due to the galaxy's proximity to the group cen- 
troid, its shear contribution at low values of fhaio is ap- 



proximately degenerate with the shear contribution from 
the group dark matter halo at large values of fhaio- The 
extreme luminosity of this galaxy may be the result of 
blending in our photometry, but since its mass is con- 
centrated over a narrow range of position angles relative 
to the lens, it is unlikely that modeling this object as 
the superposition of less massive galaxy halos will make 
a large differen c e in th e shear. 

iWavth et all (|2005l ) model Q0047 as an elliptical 
power-law model and find that the best fit is near- 
isothermal and requires no external shear to explain the 
data. It is possible that their models could accommo- 
date a small amount of shear, which would be consistent 
with our shear distribution for this system. We find that 
our shear distributions are consistent with zero shear to 
within the 95% confidence contours. 

Most of the uncertainty in the shear in MG0751 comes 
from the scatter in the FJ relation, which increases the 
error in 7 from 0.02 to 0.04. There is a very bright 
group galaxy close to the lens (I = 17.89, 0.1' separation) 
that is likel y the main contributor to this uncertainty. 
iLehar et al.l JT997) fit point mass and SIS lens models 
with external shear to this lens, and their shears are in- 
consistent with our environment shears. However, their 
models do not include both ellipticity in the lens galaxy 
and external shear, and thus are not directly comparable 
to ours. 

The host group of MG1654 has a velocity dispersion 
three times lower than most of the other groups, corre- 
sponding to a virial mass smaller by at least an order 
of magnitude and resulting in a small shear contribution 
from the group environment. The velocity dispersion of 
this system may be intrinsically low, although it could be 
underestimated due to the small number (eight) of con- 
firmed group galaxies in our spectroscopic sample and 
the lack of photometric red sequence group members. 
The small shear from the group environment also leads 
to a very small scatter in the shear due to the observa- 
tional errors we considered. As a result, most of the shear 
amplitude and scatter in our environment model comes 
fro m the LOS galax ies. A modeling analysis performed 
by lKochanekl (|1995[ ) tested a number of lens models and 
found that a quasi-isothermal ellipsoid (0.9 < a < 1.1) 
provided the best fit. The range of shears and position 
angles obtained from his models are consistent with our 
shear distributions. 

6.5. Remaining Sources of Error in Environment 
Analysis 

We have modeled the environments of our lens systems, 
accounting for various observational and theoretical un- 
certainties (§ |4.4[) . However, there are other sources of er- 
ror that could affect our results and that may explain the 
discrepancy between our shear distributions and those 
from lens modeling (§ 16. 2p . including: 

• We assume that all galaxies in our analysis are 
early-type galaxies. We apply the Faber- Jackson 
relation to obtain velocity dispersions and use an 
elliptical galaxy SED when performing bandpass 
corrections. To roughly test this assumption, we 
estimate the fraction of red galaxies in our spec- 
troscopic sample both by looking for emission lines 
in the spectra (Momcheva et al. in preparation) 
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and comparing the observed colors to those of an 
Sa-type g alaxy template fro m the GALEV2 tem- 
plate set dBicker et al.ll200"4[ ). These tests suggest 
that ~60-80%3 of galaxies in the host group of the 
lenses and ^50% of the LOS galaxies are red, so 
this assumption is not baseless. Id eally, we would 
like t o apply the Tully-Fisher (TF; iTullv fc Fisherl 
119771 ) relation to the blue galaxies, but this is diffi- 
cult because we lack information on the inclinations 
of disk galaxies. In addition, the relative errors in 
the FJ and TF relations at intermediate redshifts 
are not well-constrained enough to determine a ro- 
bust calibration between the two. While we do not 
consider evolution in the FJ relation due to veloc- 
ity dispersion evolution, we do account for passive 
luminosity evolution in the galaxies. 

iMandelbaum et all (|2006t ) find that late-type 
galaxies with luminosities L ~ to 2L* have 
roughly half the r-band mass-to- light ratio of early- 
type galaxies at similar luminosities. At fainter lu- 
minosities (L ~ L*/2), the mass-to-light ratios are 
nearly equal. As a crude check of the FJ assump- 
tion, we run a test where galaxies with blue colors 
are given half the mass (and therefore, l/y/2 times 
the velocity dispersion) that they normally would 
be assigned by the FJ relation. Blue galaxies are 
defined as galaxies with a. V — I or R — I color 
bluer than that of the Sa-type galaxy template at 
the same measured or assigned redshift. This con- 
servative test provides intuition as to how much 
our shears are affected by the assumption of the 
FJ relation for all galaxies. We find that this pro- 
cedure changes the the shear distributions by less 
than 0.01 in either component and that the mis- 
matches between our environment shears and the 
lens model shears remain. 

• In our analysis of the LOS, we take individual 
galaxy halos into account, but ignore the effects 
of group or cluster-sized halos. In principle, we 
could perform a similar fhaio analysis for each as- 
sociation along the line of sight. We test this effect 
in the six fields where we have lens model shears. 
We take galaxies in each galaxy group along the 
LOS identified by Momcheva et al. (in prepara- 
tion) and replace them with an SIS halo of mass 
determined from the calculated group velocity dis- 
persion. We cross-correlate the results of these tri- 
als with those of our standard Monte Carlo trials 
as in § 16.21 and find that the shears agree to within 
the 68% confidence contours. Accounting for the 
halos of massive LOS structures does not change 
the shear distributions appreciably, nor bring the 
shears into agreement with the lens model shears 
in the three cases where they are discrepant. 

• Our assumptions in the scaling of the group galaxy 
halo masses could be incorrect. Some galaxies near 

9 This fraction is higher than the fraction of all confirmed group 
members tha t lie on the photometr ic red sequences (Table [2] fur- 
ther details in Williams et al. ( 2006) and Williams et al. (in prepar- 
tion)), but this increase is expected given the conservative defini- 
tion (i.e., much narrower color cut) of the photometric red sequence 
galaxies f§ l4,3jl . 



the centers of groups may have had their dark mat- 
ter halos tidally stripped, while others on the out- 
skirts may st ill have extended hal os. Using N-body 
simulations, Ghigna et al. (1998]) find a marginal 
dependence of halo extent relative to virial radius 
on projected separation from the center of a host 
galaxy cluster, s ugges ting that this effect is small. 
iLimousin et al.l (|2009? ) perform a similar analysis 
and find a stronger dependence. However, is not 
clear how these analyses relate to the lower mass 
groups in our sample. 

• We assume that the groups are all virialized, but 
if they are not, the velocity dispersions are im- 
perfect measures of the group potentials. For 
two of the three systems where our results do 
not agree with the lens model shears (B1422 and 
RXJ1131), we find a photometric red sequence at 
the group redshift, suggesting a relaxed component 
(|Zabludoff fc Mulchaevlll998l) . 

• We assume all galaxies in our environment anal- 
ysis are singular isothermal spheres for simplicity, 
but in reality, they are likely to be b etter mod- 
eled a s singular isothermal e l lipsoids (iTreu et ail 
20061: iKoopmans et al l 120061: iBolton et al.l 120081: 
Koopmans et al.l [2009). This effect is probably 
small because the mass distribution of a perturber 
outside Re of the lens is still centered at the same 
projected separation and position angle from the 
lens. Changes in the shape of the perturber's mass 
distribution will only be small effects on the po- 
tential at the lens position. Furthermore, the frac- 
tional change in the shear can only be of order the 
cllipticity, and it can be an increase or decrease 
depending on how the perturber is oriented with 
respect to the lens. It would take a fairly patholog- 
ical configuration — many very elongated galaxies 
oriented in just the right fashion — to produce a 
dramatic change in the shear. 

• The group centroids determined from the galaxy 
positions could be offset from the true mass cen- 
troid. Luminosity-weighting the centroid determi- 
nation does not affect our results. However, X- 
ray gas comprises up to 10% of the total mass 
of the group for groups with velocity dispersions 
similar to the ones we consider in our sample 
(jGonzalez et al.ll2"007| ) and could be offset from the 
galaxy-determined centroid by ~ 10— lOOkpc. The 
effect of this offset on the shear will be small if the 
projected distance from the lens to the galaxy cen- 
troid is much larger than the likely offset between 
any gas component and the centroid. We run tests 
in B1422 and PG1115 where we fix the group halo 
ce ntroid at the posi t ion of the X-ray centroid given 
by iFassnacht et al.l (|2008| ) , but this results in the 
environment shears moving further from the lens 
model shears, indicating that this is not the cause 
of the discrepancy. 

• Our method for populating the LOS redshift distri- 
butions does not handle voids properly, as a given 
galaxy has an uniform probability of being assigned 
a redshift between two adjacent galaxies in a given 
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bin, and we do not have the statistics to define 
voids well. These voids could contribute a neg- 
ative shear since they are underdense relative to 
the mean density along the line of sight. Voids 
are likely less significant than overdense regions be- 
cause the deviations from the mean density along 
the LOS are smaller than prominent peaks, which 
can have arbitrarily high overdensity. Through 
Monte Carlo simulations, Momchc va et al.l ([2006) 
found that voids are unlikely to contribute enough 
negative convergence to counter that due to promi- 
nent peaks along random lines of sight. Thus, we 
infer that their contribution to the shear is even less 
significant because cancellation effects arise from 
multiple voids and shear is a tensor quantity while 
convergence is a scalar. 

6.6. Additional Sources of Error in Lens Models 

There are several potential sources of error in the shear 
values derived from lens modeling, including: 

• Modeling the lens galaxy as an ellipsoid with a sin- 
gle power-law density profile may not be accurate. 
Two-component model fits to the density profile, 
encompassing both a stellar component and dark 
matter compon ent, may be more rea listic in cer- 
tain cases (e.g. iKochanek e t al. 2006]), as they al- 
low for a more complex angular structure. We fit a 
de Vaucouleurs + unconstrained NFW model (rep- 
resenting the stellar and dark matter components, 
respectively) to B1422 and PG1115 to test whether 
the shears are significantly different. While the 
shear distributions shift by ~ 0.1 compared with 
the power-law models, they result in PG1115 be- 
coming inconsistent with our environment shears 
and do not result in agreement between the model 
and environment shears in B1422. 

• The lens galaxy may be tidally truncated, which 
would change its assumed mass distribution. The 
effect of this truncation is probably small because 
Re for lens galaxies is typically ~10 kpc, which 
is consistent with an effective radius for a massive 
galaxy. 

• Substructure lensing, including substructure in the 
lens galaxy, the group environment, and along the 
line of sight that is projected within Re of the 
lens galaxy, is not accounted for in the lens mod- 
els. This substructure co uld introduce perturba- 
tions in the lens potential dMao fc S chneider 1998; 
IMetcalf fc Madaull200lt IDalal fc Kochanekl I2002T ) . 
resulting in astro metric perturbatio ns in the lensed 
image positions ()Chen et al.l 120071) . Preliminary 
work on HE0435 suggests that substructure does 
not significantly shift the model shear distributions 
(Fadely et al. in preparation), but all the systems 
need to be examined in detail to rule out substruc- 
ture as a significant systematic effect in lens model 
shears. 

• The lens models do not account for higher order 
terms than 7 in the expansion of the lens poten- 
tial. Higher-order perturbations may result in bet- 
ter fits to the data with different values of 7 than 
are currently being calculated. 



• Like most previous lens modeling analyses, our lens 
models omit non-linear couplings between the lens 
plane and line-of-sight perturbers. Including those 
non-linear effects might yield better fits with differ- 
ent values of 7. It might also have the consequence 
of modifying the "flavor" of shear that provides the 
best comparison between lens models and the envi- 
ronment analysis (although, as noted in § 14.21 the 
choice of shear flavor does not significantly affect 
our environment-derived shear distributions). We 
are currently exploring the impact of different shear 
flavors on lens models (Keeton et al. in prepara- 
tion). 

• The lens models assume k = 0, as standard anal- 
yses do. This assumption is incorrect, as k > in 
these systems (Wong et al. in preparation). We 
test the effects of this assumption by calculating 
the reduced shear, g c s = 7 C , S /(1 — k), from the 
7 and k that we obtain from the lens environment 
analysis and comparing the two reduced shear com- 
ponents to the shears from the n = lens model. 
Our results do not change: the mismatches remain 
discrepant at the > 95% level and the matches are 
the same as before. As a result, and because we 
wish to define the environment-determined shear 
as the true shear, independent of the convergence, 
we tabulate and plot the true shears as our shears 
throughout this paper. 

7. CONCLUSIONS 

Using new spectroscopic and photometric data from a 
survey of nine gravitational lenses in groups of galaxies, 
we quantify the effect of the local and line-of-sight envi- 
ronments on the lens potentials by directly constraining 
the environmental shear. We also analyze the relative 
importance to the shear of observational and theoretical 
uncertai nties, including uncertai nty in the Faber- Jackson 
relation (|Faber fc JacksonH l976). errors in the projected 
group centroid position and velocity dispersion, galaxy 
magnitude errors in the photometric data, uncertainties 
in the form of the mass profile of the group halos, and 
uncertainty in the concentration parameter of the group 
halos. The full lens environment contributes a significant 
shear of 7 = 0.08 on average, ranging from 0.02 to 0.17. 

For the six 4-image lenses where we can compare our 
shears to those derived from lens modeling, our envi- 
ronment analysis can reproduce the lens model shears 
in HE0435, PG1115, and HST14113. However, for the 
other three systems (RXJ1131, B1422, WFI2033), the 
environment shears are inconsistent with the lens model- 
derived values at more than the 95% level, pointing to a 
serious problem in the lens modeling or in the way that 
we have characterized the environment. There do not 
appear to be any characteristics common among the sys- 
tems where our shears are inconsistent with lens models, 
nor do the models themselves necessarily have large x 2 
values relative to the number of degrees of freedom. 

The contribution to the shear from structures along 
the line of sight to the lens can have an effect on the 
order of that produced by the local group environment. 
The average environmental shear is 7 = 0.08, compared 
to an average shear of 7 = 0.05 when we only consider 
perturbations from the local environments of the lenses. 
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The mean offset in (j c , 7,5) space when the effects of LOS 
perturbers are added is Aj env = 0.06. Most of the LOS 
shear comes from objects projected within ~2' of the lens 
and brighter than / = 21.5. Less than 0.01 in either 7 C 
or 73 is contributed on average by galaxies > 5' away 
or fainter than / = 21.5, so these generally can be ne- 
glected for lenses at these redshifts. We conclude that to 
minimize uncertainties in the shear from spectroscopic 
incompleteness, future spectroscopic surveys of the en- 
vironments of lenses at these redshifts should prioritize 
/ < 21.5 objects projected close (within ~2') to the lens. 

Our analysis shows that accounting for many possi- 
ble observational and theoretical uncertainties typically 
shifts the shear by ~0.01 in either 7 C or 7 S . However, 
these uncertainties add ^0.015 to the scatter in our shear 
distributions on average. If these uncertainties are not 
taken into account, the scatter in the individual shear 
components ranges from 0.01 to 0.04 due to systematic 
effects alone. 

Individually, scatter in the Faber- Jackson relation (FJ; 
iFaber fc Jackson! 1 19761) and the error in the group cen- 
troid position contribute most to the scatter in our 
shear distributions, adding 0.03 to the scatter in ei- 
ther 7 C or 7 S in the most extreme cases. Improving 
the errors introduced by the FJ relation requires a scal- 
ing rel ation with less scatter such as the Fundamental 
Plane (iDiorgovski fc Davis] Il987t iBernaxdi et aT1 l2003ri 
ICappellari et al.l 120061 : iRobcrts on et al.l 120061 ). Obtain- 
ing effective radii for the galaxies in our sample would 
give better estimates of t he galaxies' internal y elocity 
dispersions. For example, iBernardi et al.l (|2003b| ) claim 
an error of 0.05 in log a and 0.01 in log Jo for the Funda- 
mental Plane, compared to the intrinsic scatter of ~0.07 
in log a that we assume for the F J relation. For a trun- 
cated SIS perturber, 7 cx a 2 inside r t and 7 cx a 3 outside 
r t , so a reduction in the scatter in a can have a large 
effect on the inferred shears. 

The errors in the group centroid position are most sig- 
nificant for systems where the lens is close to the centroid, 
particularly HE0435, where one of the shear components 
changes by 0.05. More complete spectroscopic sampling 
of group members and X-ray imaging of the group's hot 
gas distribution could reduce the centroid errors. 

Magnitude errors in our photometric data contribute < 
0.01 to the scatter in either shear component. The errors 
in the group velocity dispersion propagate into errors in 
the total mass of the host groups, typically adding ^0.01 
to the shear scatter. The errors in the concentration 
parameter of the group dark matter halo are similarly 
small. 

We do not know how the mass in a group is appor- 
tioned between the group dark matter halo and the in- 
dividual group galaxies because it is dependent on the 
degree to which the group members have been tidally 
stripped via galaxy-galaxy interactions as the group 
evolves. Therefore, we leave fhalo, the fraction of the 
group's virial mass in the common group halo, as a free 
parameter. Yet, even varying fhalo does not allow us to 
reproduce the lens model-derived shears in half of our 



4-image lens subsample. 

Another theoretical uncertainty is the mass density 
profile of the group dark matter halo. To test the effects 
of this uncertainty, we have analyzed the shear profiles 
of a SIS and NFW group halo. The difference in the 
shear between the two mass profiles is at most ^0.02 ex- 
cept when the lens is projected within ~ r s /2 of the halo 
centroid, assuming large values of fhalo- This difference 
between the mass profiles may have a small but notice- 
able effect for lenses with a projected separation from 
the group centroid near this peak, but the choice of halo 
profile appears to be otherwise unimportant. 

The disagreement between the shears calculated di- 
rectly from our environment analysis and those derived 
from lens modeling might arise from the fact that lens 
modeling generally ignores higher-order expan sions of 
the lens potential (jKeeton fc Za bludoff 120041 ) . makes 
simplistic assumptions about the form of the lens galaxy 
mass pr ofile, and does not account for substructure tens- 
ing (e.g.lMao fc Schneide ^ll99llMetc alf fc Madaul l2Q0H 
Dal ai fc Kochanekl I2002T There are also possible bi- 
ases or sources of error that we have not considered in 
our environment analysis. These effects include location- 
dependent tidal stripping of galaxy halos and variations 
in galaxy morphology. We have now quantified the envi- 
ronmental perturbations on strong gravitational lens po- 
tentials further than previous studies. Future work must 
characterize these other sources of uncertainties to fully 
understand the limits of gravitational lens constraints on 
important quantities like Hq and the properties of lens 
galaxy halos. 

Our methodology also permits us to calculate the con- 
vergence (k) introduced by the environment (Wong et 
al. in preparation). Unlike the shear, k is not con- 
strained by the lens models. While statistical constraints 
on k may be inferred from methods s uch as ray-tracin g 
through cosmological simulations (e.g. lSuvu et alfe OlO). 
measurements of the lens environments are the only way 
to directly determine its effect on lens-derived quanti- 
ties. Were the convergence values for lens systems com- 
parable to their shears (i.e. ~0.1), Hq determinations, 
whos e errors scale as (1 — k), w ould be significantly bi- 
ased (Kceton & Zabludoff 2004). Our work shows that 
lcnsing constraints must consider both local environment 
and line-of-sight galaxies in their error budgets. 
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APPENDIX 

A. SHEAR CALCULATION FORMALISM 

We present the details of o ur for malism for calculating the shear due to the host group halo (Appendix IA. II) and 
individual ga laxies (Append ix IA.2|) . The calculations are performed using an updated version of the lensing software 
developed by lKeeto n (2001). For the code to perform the necessary shear calculations, we must input parameters that 
we calculate based on our data and various assumptions that we detail in this section. 

A.l. Shear due to SIS and NFW Group Halos 

We test how large an effect the choice of an SIS or NFW profile for the group halo has on the lens potential. The 
shear due to an SIS halo is 

is is = ^ , (Al) 

where R is the angular offset of the halo centr oid from the lens, while Re is the Einstein r adius given by Equation [§1 
To calculate the shear due to an NFW halo ([Bartelmannl 119961 : IWright fc BraineTdl 120001) , we need to calculate the 
parameters r s and n s . The NFW scale radius r s is that at which the power law slope of the density profile is equal to 
that of an isothermal sphere. r s is often expressed in terms of the concentration parameter, c V i r — r V i r /r s , where r V i r 
is the virial radius of the halo. The dimensionless parameter n s is defined as 

«. = ^ (A2) 

where p s is the central density of the halo and E c is the critical surface density for lensing, 

c^Ds 
c AttGD ls D l ' 

Dl, Ds, and Dls are the angular diameter distances between the observer and the lens, the observer and the source, 
and the lens and the source, respectively. Using the radial velocities of the spectroscopically confirmed group galaxies, 
we estimate the velocity dispersion, a grp , and the mass of the dark matter halo, M/> a ; , of each group. In the shear 
analysis, we add photometric red sequence galaxies to the group, and these quantities are redetermined for each 
Monte Carlo trial using the bootstrap method as noted in the text. For this test, however, we assume that only the 
spectroscopically confirmed members are in the group. We then calculate the parameters needed for the shear analysis. 
We first approximate the virial mass of the group, M V i r as 

M mr = -nrl ir A c (z)p c (z), (A4) 

where p c {z) and A c (z) are the critical density and the characteristic overdensity at the lens redshift ([Eke et al.lll998[) 
respectively. From our assumed cosmology, we have 

_ ^jmu^f+M (A5) 

07TG- 

The characteristic overdensity, A c , is the ratio of the mean density inside t he group halo to the cr itical density at the 
redshift of the group. We use a form of the characteristic overdensity from Eke ct al. (1998, 2001): 



(A3) 



A r = 178 



= 178 



/ nt 0.45 
Pm(Z) 



Pc(z) 



n m (i + zr + n A \ • (A6) 

This form is accurate to within 5% of a spherical collapse model ([Eke et al.lll998D for our cosmology. We account for 
this uncertainty by allowing for a 5% Gaussian scatter in the characteristic overdensity, which propagates directly into 
a 5% error in the virial mass. We test the effects of this uncertainty explicitly in our shear analysis and find that it is 

negligible in comparison with the other sources of error. 

We assume the virial theorem and express the virial mass (e.g. iKoranvi fc G cllcr 20Q3) as 

Mmr = (A7) 



We then equate this to M v i r from Equation IA4I and solve for r V i r as a function of a 



'grpi 

3ggrp (A8) 



2^/irGA c p c 



Now we are able to evaluate r V i r in terms of known or assumed quantities and can calculate Mhaio for an assumed 
NFW halo from Equation IA41 To determine the concentration parameter c V i r for the group halo, we use the "ENS" 
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code described in lEke et all (|2001l ). When running this code, we assume erg = 0.8 and a shape parameter T = 
(n^e-^-V^^"/^) ra 017 in our cosmology (|Sugivamalll995h . 

We now determine r s for each group by setting r s — r V i r /c V i r . To determine p s , we set Mhalo equal to the NFW 
density profile integrated over a spherical volume enclosed by the virial radius. Solving for p s , we find 



_ Mhalo 

Ps ~ 4irr s 



-i -1 



(l+r/r s y 

„3 



47rr^ ir [ln(l + c mr ) - c mr j (1 + c mr )] 



(A9) 



and calculate k s . Using these parameters, we determine the shear profiles for the SIS and NFW halos for each group. 
This calculation assumes that all of the group mass is in the dark matter halo, but we can rescale this by changing 
Mhalo to some fraction of the total virial mass as in § 15.1.31 This is not exactly correct since the concentration of the 
halo determined from simulations implicitly accounts for the individual galaxies as subhalos. However, the fractional 
change in the concentration between trials where Mhn.io = M v ^ and those where M hnin << M V i r is generally smaller 
than our assumed scatter of log c V i r = 0.14 (jBullock et al.ll200ll :IWechslcr et al. 2002), so the rescaling of the halo mass 
gives reasonable results. 

Figure [7] shows the dependence of shear on projected offset from the halo center of mass for each of our nine groups. 
The SIS and NFW halos are assumed to contain all the group mass, which has been calculated from the mean velocity 
dispersion determined across multiple Monte Carlo trials. The absolute value of the difference between the SIS and 
NFW shear profiles (| Ajhalo] m the inset panels in Figure [7]) show that the choice of profile does not make a significant 
difference except where the lens is projected within < r s /2. There is a small local peak in |A7^ ; | at larger radius, 
but its magnitude is ^0.02 or less. For low (< 0.5) values of fhaio, the typical shear difference at this peak would be 
roughly ~ 0.01 or less. These results are consistent with the results of iWright fc Brainer d (2000), who find that the 
mean shear ratio between NFW and SIS halos within a virial radius for a ACDM cosmology varies by w25% for halo 
masses like ours. 

A. 2. Shear due to Truncated Singular Isothermal Sphere Galaxy Halos 

In our analysis, we model the mass distribution of a galaxy-scale perturber as a truncated singular isothermal sphere 
(TSIS). For a general SIS with velocity dispersion a, the shear amplitude is given by Equation I AIL When working 
with TSIS profiles, we use a modified version of the basic SIS formalism. The density profile of an isothermal sphere 
with velocity dispersion a is 

PsMr) = ^ (A10) 

Suppose we want to truncate this profile at some radius r*. One way to obtain a smooth truncation is to write the 
density profile as 



for some n > 0. Larger values of n correspond to sharper truncation. We can compute the lensing properties of a 
TSIS perturber for integer values of n. For n — 6, we find the convergence k, shear amplitude 7, and deflection angle 
a to be: 

= If ~ ff^t 25 ^ 1 + 128r3 + 96r5 + 8 ° r7 + 7 ° r9 + 63 £~ U ]' (A12) 

i(R) = |f - |f^[i26(i + er 1 + isor 1 + 2r 3 - 3or 5 - ^r 7 - ^c 9 - ear 11 ], (ais) 

a{R) =R E + h256£ + 63 + 128r 1 + 32£~ 3 + 16r 5 + 10C 7 + ?C 9 ] - ( A 14) 

2bbKl)p 



where £ = v/l + R 2 D 2 p /rf and Dp is the angular diameter distance to the perturber. We assume n = 6 throughout 
our analysis when discussing TSIS profiles since it represents a sufficiently sharp truncation and larger values do not 
make a noticeable difference in our results. 

B. TRUNCATION RADII OF GROUP GALAXIES 

When apportioning mass among the group halo and the group galaxies, the group galaxies' truncation radii are 
scaled so that the density at the truncation radius is the same for all group galaxies. If we assume infinitely sharp 
truncation, a sing ular i sothermal sphere with a truncation radius rt and velocity dispersion a has a density profile 
given by Equation IA10I out to r t . Beyond r t , there is no mass, so p(r > r t ) — 0. The galaxy has a total mass of 

Msis=—q 1 - (Bl) 
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Fig. 7. — Shear profiles for SIS (solid line) and NFW (dashed line) group halos. The halos are assumed to contain all of the group mass. 
We plot the NFW scale radius r a (dotted line) and the projected offset of the lens galaxy from the group centroid (dash-dotted line). The 
inset panels show the absolute value of the difference in the shear profile between SIS and NFW group halos. The peak shear difference 
between the two profiles is roughly ~0.02 except for r < r s /2. Only HE0435 and PG1115 are close enough to the group centroid to fall 
within this region. 



Given a group with a virial mass M V i r , a halo mass fraction fhaio (Equation 1 10[) . and iV galaxies with masses Mi and 
internal velocity dispersions cr^, we want to find the truncation radii r t i such that p{r t i) is the same for all galaxies. 
We do this by setting the total mass in the group galaxies equal to the sum of the galaxy masses, 

N 

(1 - fhalo)M mr = Mi 
i=l 

2 N 
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where p t is the density at the galaxies' truncation radii. Solving for this density, we hnd that 



-i 2 



2 



E^ 3 



(B3) 



Pt = 



*G 3 0— fhalo)*MZ 



vir 



i=l 



and the corresponding truncation radius of a given galaxy j is 

.. . - a i 



-, -l 



= — fhalo)M v ir ^ of 



(B4) 



i=i 



C. EFFECTS OF RADIAL AND LUMINOSITY SAMPLING ON LINE-OF-SIGHT SHEAR 

C.l. Projected Separation Cut 

Our spectroscopic catalog for each system includes several hundred objects projected within ~15' of the lens. How- 
ever, our photometric catalog includes objects within ^25' of the lens and goes to fainter magnitudes, containing ^10 4 
objects per system. With this many objects, it is not possible to achieve reasonable computation times for our Monte 
Carlo simulations, where we generate and analyze thousands of redshift realizations of the entire ensemble of galaxies 
in our catalog. However, most of these galaxies are likely to have a negligible effect on the shear. Therefore, we cut 
all LOS objects with a projected separation >5' from the lens, reducing the number of galaxies in our photometric 
catalogs by over an order of magnitude. 

We expect that the shear due to LOS perturbers will decrease with projected separation from the lens for several 
reasons. First, the shear from an individual perturber with an SIS profile (or any other reasonable density profile that 
drops off with radius) will decrease with projected separation. Given that shear adds as a rank-2 tensor quantity, we 
also expect that the total shear amplitude due to all objects interior to some radius from the lens will approach some 
limiting value as we include more perturbers (assuming that their azimuthal positions around the lens are random). 
This is the result of canceling effects from perturbers lying on lines orthogonal to each other in sky coordinates centered 
on the lens. Within larger radii, we are enclosing more LOS objects, so the shear contribution from all objects near 
the edge of the reg ion will tend to be smaller than objects closer to the lens where shears are less likely to cancel out. 
lAuger et al.1 (j2007j ) investigated environment effects on two lenses in galaxy groups and found that their contributions 
were dominated by galaxies projected within 15" of the lens, which supports our justification for a cut on projected 
separation. 

Despite these effects, one could imagine an unlikely configuration of LOS objects arranged as to create a large shear. 
To test the validity of a projected separation cut, we run our redshift randomization procedure on objects in bins of 
projected se parat ion from the lens. We perform this test for objects brighter than I — 21.5, a limit that we justify 
in Appendix IC.2I For each lens, we perform 1000 trials for each of six circular apertures extending 0' to 8' from the 
lens in projected separation, each with a cutoff radius 1' larger than the previous aperture. In each case, we account 
for all objects in the local group environment, but exclude LOS objects with a projected separation outside the cutoff 
radius. We then calculate the mean 7 C and j s from the objects within each aperture (Figure [8]). Based on the amount 
of variation with cutoff radius, we conclude that most of the LOS shear comes from objects projected within 2-3' of 
the lens. The mean -f c and 7 S contribution from the ob jects beyond 5' a cross all nine lens environments is smaller than 
0.01, well within the error bars. This is consistent with iBrainerdl (|2010f ) . who finds that the shear due to galaxy-galaxy 
weak lensing extrapolates to zero at a separation of ~5'. Based on this analysis, we assume that 5' is a reasonable, 
conservative cut-off in projected separation. 



For LOS objects in our photometric sample that are fainter than our spectroscopic limit of I = 21.5, our method 
for creating a redshift distribution based on objects in our spectroscopic catalog (§ I4.3[) breaks down because we do 
not have a representative spectroscopic sample from which to draw redsh ifts. Objects in our p hotometric catalog can 
be as faint as / ~ 24, but the star-galaxy separation routine described in lWilliams et all (|2006l ) starts to break down 
around / ~ 21.5. One solution, which we employ, is to cut objects fainter than / = 21.5 from our LOS sample. The 
main assumption is that fainter objects contribute little to the shear in comparison to brighter objects. Objects fainter 
than this limit are either less massive than the typical galaxy at the low-mass end of our spectroscopic sample, or at 
higher redshift than the lens. Both effects will reduce their shear contribution. 

To justify the I — 21.5 cut, we plot the shears from the local group environment plus all LOS galaxies within 5' of 
each lens as a function of cutoff magnitude (Figure |9]). In each bin, we only include the shear from LOS galaxies 
brighter than some limiting magnitude all the way down to / = 22.5. For the bins fainter than / = 21.5, we draw 
redshifts from our faintest galaxy redshift distributions (roughly I = 21 - 21.5) in § 14.31 We run 1000 trials for each bin, 
plotting the mean and standard deviations of j c and 7 S . Because our redshift assignment procedure draws redshifts 
based on a distribution of galaxies limited at I = 21.5, it will tend to underestimate the masses of galaxies fainter 
than this limit. However, this effect will be offset to some extent by the fact that their true redshift distribution is 
likely shifted away from the lens redshift to higher redshifts, reducing their shear contribution. 
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Fig. 8. — Shear contribution from the entire local group environment plus only those LOS objects brighter than / = 21.5 and within the 
given projected separation from the lens. Very few objects in the local environment are fainter than / = 21.5 or lie outside 5'. The points 
represent the mean shear and standard deviation over 1000 trials for 7 C (blue) and -y s (red). The 7 S points are offset by 0.1' for clarity. 
The first bin represents the shear from the local group environment alone. Our projected separation cutoff is 5' (dashed line). Most of the 
variation in the shear comes from objects projected within 2' of the lens, whereas objects projected further than 5' away contribute less 
than 0.01 to the shear on average. 
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Fig. 9. — Shear contribution from the enitre local group environment plus only those LOS objects projected within 5' of the lens and 
brighter than the given / magnitude. Very few objects in the local environment are fainter than / = 21.5 or lie outside 5'. The points 
represent the mean shear and standard deviation over 1000 trials for f c (blue) and ■js (red). The 7 S points are offset by 0.1 mags for clarity. 
The first bin represents the shear from the local group environment alone. Our magnitude cutoff is / = 21.5 (dashed line). Most of the 
variation in the shear due to the LOS comes from objects brighter than 21.5, while fainter objects contribute less than 0.01 to the shear on 
average. 
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The results of Figure [S] show that much of the LOS shear comes from bins brighter than I = 21.5 (our spectroscopic 
limit). At fainter magnitudes, the difference in the shear components is ~0.01. As with the projected separation, 
there tends to be a decline in additional shear contributions as we go to fainter magnitudes, and the variations are well 
within our error bars beyond I = 21.5. For this analysis, we adopt this conservative magnitude cut to be consistent 
with the limit of our spectroscopic catalog. 
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